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ABSTRACT 

We present key sample results of a systematic survey of the pulsation properties of models 
of hot B subdwarfs. We use equilibrium structures taken from detailed evolutionary sequences 
of solar metallicity (Z = 0.02) supplemented by grids of static envelope models of various 
metallicities (Z = 0.02, 0.04, 0.06, 0.08, and 0.10). We consider all pulsation modes with 
1 = 0, 1,2, and 3 in the 80-1500 s period window, the interval currently most suitable for 
fast photometric detection techniques. We establish that significant driving is often present 
in hot B subdwarfs and is due to an opacity bump associated with heavy element ionization. 
We find that models with Z > 0.04 show low radial order unstable modes; both radial and 
nonradial (p, f, and g) pulsations are excited. The unstable models have T eS < 30,000 K, and 
log g < 5.7, depending somewhat on the metallicity. We emphasize that metal enrichment 
needs only occur locally in the driving region. On this basis, combined with the accepted view 
that local enrichments and depletions of metals are common place in the envelopes of hot B 
subdwarfs, we predict that some of these stars should show luminosity variations resulting 
from pulsational instabilities. 

Subject headings: stars: interiors— stars: oscillations— subdwarfs 
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1. INTRODUCTION 

In the last twenty years, considerable progress has been 
made in our understanding of the physical properties and 
evolutionary status of hot, hydrogen-rich subdwarf B stars 



(see, e.g., preenstein & Sargent 1974 


; Heber 1987 




1991; 


Allard et al. 1994 




Saffer et al. 1994; Villeneuve et 



al. 1995). Following the sce nario originally pro posed by 
Heber and collaborators (e.g., Heber et al. 1984 ), it is cur- 
rently believed that subdwarf B stars are ~ 0.5 Mq objects 
with hydrogen envelope masses that are small enough that, 
after core helium exhaustion on the extended horizontal 
branch (EHB), they never evolve toward the asymptotic 
giant branch (AGB). They remain at high effective tem- 
peratures (T c ff >, 20,000 K) throughout their core-burning 
evolution. In the post-HB, He-shell-burning stage of evo- 
lution the models are referred to as 'AGB-Manque" ob- 
jects (see Greggio & Renzini 1990 ) and are associated 
with the field subdwarf O stars (Dorman, O'Connell & 
Rood 1995 and references therein). A number of post- 
EHB evolutionary studies have found their way in the lit- 
erature in the last several years (see porman 1995 for a 
review). On the basis of the evolutionary models of Dor- 
man, Rood, & O'Connell 1993, |Bergeron et al. 1994| have 
shown that the subdwarf B stars most likely evolve to be- 
come the (relatively rare) low-mass DAO white dwarfs, 
thus contributing to a small percentage of the total white 
dwarf population. 

While asteroseismology is proving to be an extremely 
powerful tool in studying other types of stars, its poten- 
tial has not yet been studied for subdwarf B stars partly 
because of the lack (until recently) of sufficiently realistic 
equilibrium structures. Furthermore, and to our knowl- 
edge, luminosity variations caused by pulsational instabil- 
ities have not been reported for these stars^j. However, pul- 
sation theorists with a keen eye may have noticed for some 
time that the potential for driving pulsation modes appears 
to exist in the envelopes of s ubdwarf B stars. Indeed, early 
models of such enve lopes (Winget & Cabot 1981 ; We- 
semael et al. 1982; IVilleneuve 19821; Groth, Kudritzki, 



& Heber 1985) as well as their more sophisticated mod- 
ern counterparts are all characterized by the presence of a 
Hell-Helll convection zone. By analogy with other types 
of pulsating stars (whose instabilities are always driven by 
one form or another of an opacity mechanism), one might 



2 Note that, in th e course of their Stromgren photometric program. B erg- 
eron et al. 1984 noted and subsequently reported the possible variabil- 
ity of the classical subdwarf B star Feige 108. However, the equipment 
available to them as well as time constraints did not allow the gathering 
of a light curve of sufficient quality for time-series analysis. 



expect that the opacity bump associated with this partial 
ionization zone could also excite pulsation modes in hot B 
subdwarfs. 

Motivated in part by this observation (and see Van Horn 
1991), we undertook a systematic exploration of the po- 
tential of asteroseismology for subdwarf B stars. The pur- 
poses of the present Letter are (1) to report on the salient 
results of our investigation, (2) to call attention to the very 
real possibility that some subdwarf B stars may undergo 
stellar oscillations, and (3) to encourage searches for lu- 
minosity variations in such stars. 

2. COMPUTATIONS AND RESULTS 

The first batch of equilibrium models investigated in 
this study consists of full stellar models taken from 5 dis- 
tinct evolutionary sequences. These were chosen in or- 
der to map a significant fraction of the region of the Teg- 
log g plane actually occupied by the known subdwarf B 
stars (see Saffer et al. 1994 and references therein). In 



this particular region, we retained about 25 models per se- 
quence. The evolutionary models were computed with the 
methods described in Dorman (1992a,b) and Dorman et al. 
19 93. The new models use the OPAL opacities described 
by Rogers & Iglesias 1992 computed in December 1993, 
which adopted the element mix referred to as 'Grevesse 
& Noels 1993'. Where necessary (during He-flashes), we 
used new low temperature opacities by D. R. Alexander 
(1995, private communication, described in Alexander & 
Ferguson 1994) which were computed with the same el- 
ement mix. These smoothly match the OPAL opacity set 
within the hydrogen ionization zone. The other difference 
in the physics was the use of the Itoh et al. (1983; 1993a,b; 
1994a,b) conductive opacities. 

Each sequence describes the evolution of an AGB- 
Manque star with a core mass of 0.4758 Mq . The se- 
quences differ in that different envelope masses are con- 
sidered: Menv = 0.0002, 0.0012, 0.0022, 0.0032, and 
0.0042 Mq. The composition of the envelopes is assumed 
solar (X = 0.70388, Z = 0.01718). The luminosity of 
each model is provided by He burning confined to a very 
small core. Residual H burning at the base of the H-rich 
envelope contributes negligibly to the luminosity until the 
AGB-Manque phase of the evolution. 

We first carried out a detailed adiabatic survey of the 
evolutionary models with the help of the Galerkin finite- 
element code of Brassard et al. 1992. Specifically, for 



each equilibrium model considered, we computed pulsa- 
tion periods in the adiabatic approximation for all modes 
with 1 = 0, 1, 2, and 3 in the 80-1500 s period win- 
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dow. This interval corresponds to the range of periods 
most easily detectable with present-day fast photometric 
techniques. Moreover, it is well known that surface can- 
cellation effects render detection of luminosity variations 
excessively difficult for I > 4 in compact stars (Dziem- 
bowski 19 77; Robinson, Kepler, & Nather 198$ Bras- 
sard 1987; Brassard, Fontaine, & Wesmael 1995), so such 



modes are not of direct observational interest. We will re- 
port elsewhere the detailed results of this extensive survey, 
including discussions of the effects of changing model pa- 
rameters, the period evolution, the rates of period changes, 
and the effects of mode trapping and mode confinement 
caused by the composition interfaces in our doubly strati- 
fied equilibrium models. For the needs of the present pa- 
per, in order to give an idea of the modes to be expected in 
the 80-1500 s period window, we simply report on some 
sample results. 

Table || lists the pulsation periods for a typical model 
of a hot B subdwarf. This model belongs to the sequence 
with M coro = 0.4758 M and M cnv = 0.0012 M Q . It has 
an age of ~ 8.4 x 10 7 yrs (time elapsed since the zero age 
HB phase), a surface gravity log g = 5.46, and an effective 
temperature T e g = 27,500 K. The table gives the pulsation 
period (expressed in seconds) as a function of the radial 
order k. We distinguish between the radial modes (I = 0) 
and the p and g branches for nonradial modes with 1=1, 
2, and 3. Our results indicate that subdwarf B stars have 
rich period spectra, easily accessible with current-day ob- 
servational techniques (white light fast photometry, for ex- 
ample). Provided we can demonstrate that some of these 
modes can be excited, this should motivate observational 
searches for pulsational instabilities in these stars. 

In the second step of our investigation, we therefore 
carried out a stability analysis of our equilibrium models 
with the help of the finite-element nonadiabatic pulsation 



code briefly described in Fontaine et al. (1994) and Bras 



sard, Fontaine, & Bergeron 1996. In order to understand 
the results and be able to identify the regions of driving 
(if any), we considered the variations with depth of the 
thermal timescale and of the so-called work integral. Fig- 
ure 1 illustrates a typical case and refers to the equilib- 
rium model whose periods are provided in Table ^ The 
solid curve in the upper panel shows how the local thermal 
timescale, r t h (expressed in seconds), varies as a function 
of the fractional mass depth, logg = log(l — M(r)/M), 
in this particular unperturbed model. As a comparison, 
the values of the periods of the p(k — 3), f(k = 0), and 
g(k = 3) modes with I — 2 are indicated. In addition, the 
locations of the He burning core, the H burning shell, and 
the Hell-Helll convection zone are also shown. Because 



the local thermal timescale in the thermonuclear burning 
regions is orders of magnitude larger than the periods of 
interest (in addition to the fact that the H burning shell is 
extremely weak), it is clear that those regions cannot drive 
efficiently the pulsation modes. In practice, the e mech- 
anism thus appears to be irrelevant for exciting pulsation 
modes in hot B subdwarfs. In contrast, the local thermal 
timescale is of the same order of magnitude as the inter- 
esting pulsation periods in the envelope convection zone, 
which suggests, as anticipated above, that driving could be 
efficient in that zone. 

A look at the lower panel in the figure reveals, how- 
ever, a small surprise and an important result. The solid 
curve refers to the integrand of the work integral for the 
I = 2, / mode mentioned above. (We emphasize that, 
while specific to this particular mode, the results shown 
here are also typical of all other modes of interest.) The 
values of dW/ dr are positive in driving regions and nega- 
tive in damping regions. We find that driving is negligible 
in the convection zone. This constitutes the small surprise. 
In retrospect, it appears that the convection zone is located 
too high in the envelope, in a region that contains very 
little mass and, therefore, carries practically no weight in 
terms of driving (or damping). The important result is that 
the main driving region in hot B subdwarfs is associated 
with an opacity bump due to heavy element ionization. 
This is shown by the dashed curve which illustrates the 
run of the Rosseland opacity in the unperturbed model. 
As expected, the opacity shows a maximum in the convec- 
tion zone due to Hell-Helll partial ionization, but it also 
shows a secondary maximum around logg ~ —9.2, obvi- 
ously associated with the driving region, and which tends 
to disappear when the heavy element content (assumed to 
be Z ~ 0.02 in this model) is decreased. While the local 
thermal timescale in the driving region is larger than the 
periods of interest (see upper panel), our results indicate 
that it is not so much larger as to render the coupling to- 
tally inefficient. Moreover, the eigenfunction amplitudes 
are relatively large in the driving region. The realization 
that driving is related to the metallicity provides an impor- 
tant clue as to the possible existence of real subdwarf B 
pulsators. 

The specific mode considered in the lower panel of Fig- 
ure 1 turns out to be stable; the damping region below the 
driving region contributing somewhat more to the overall 
work integral. This result is quite typical of all of the evo- 
lutionary models computed with Z ~ 0.02. Indeed, in our 
detailed nonadiabatic survey of those equilibrium models, 
we found no unstable mode. However, in many cases, we 
found that the overall damping rates were very small, sug- 
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Table 1 

Pulsation Periods for a Typical Subdwarf B Star Model 
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gesting to us that such modes could perhaps still be excited 
in real stars (whose detailed structure might, of course, 
be different from that of the particular models we consid- 
ered). In view of the correlation we uncovered between 
driving and metallicity, and the fact that local variations of 
metallicity are expected in subdwarf B stars (see below), 
it seemed natural to investigate the effects of changing the 
metallicity. 

In the final step of our analysis, we thus constructed 
a second batch of equilibrium models taking into account 
variations of Z in the H-rich envelope. Contrary to our 
previous models, those are static (i.e., non-evolving) struc- 
tures. The models were built with a version of the stellar 



code described by Brassard & Fontaine (1994^ and Bras 



sard et al. (1996) suitably modified to produce envelope 
structures extending as deep as logq = —0.05. The con- 
stitutive physics used is nearly identical to that used previ- 
ously in the construction of our full evolutionary models, 
apart from the use of slightly newer OPAL opacities com- 
puted in 1995. Tests indicate that our ignoring of the small 
He burning core increases somewhat the values of the pe- 
riods of the g-modes (by 10-20 %, typically), but does not 
affect those of the p-modes which, contrary to the former, 
are formed essentially in the envelope region. We will pro- 
vide elsewhere more details on these models. 



We constructed 5 different model grids, one each for a 
fixed metal content Z = 0.02, 0.04, 0.06, 0.08, and 0.10 in 
the H-rich envelope (X = 0.70). Each grid consists of 72 
unperturbed models covering a rectangular region in the 
T e ff-log g plane overlapping with the region where hot B 
subdwarfs are actually found (4.34 < logT e ff < 4.62 in 
steps of 0.04, and 4.8 < \ogg < 6.4 in steps of 0.2). 
In these exploratory calculations, we fixed the total mass 
of each model to 0.48 M Q and the H-rich envelope mass 
to logq(H) = —4. Tests indicate that stability does not 
depend on the value of log q(H ). 

In agreement with our expectations concerning the ef- 
fects of metallicity, we found unstable pulsation modes for 
models with Z > 0.04. For these models, there is a blue 
edge temperature located somewhere between logT c ff = 
4.46 and 4.50, and which shows a weak sensitivity to 
metallicity. We tentatively interpret this as the conse- 
quence of the fact that metallicity affects primarily the 
size of the opacity bump (and, hence, the magnitude of the 
driving), but not strongly its location in a star of a given 
effective temperature. The blue edge itself must be due 
to the outward displacement of the driving region to the 
outermost layers (which carry little weight in the overall 
stability of a mode) as the effective temperature increases. 
Note that our actual value of the blue edge temperature 
(T c ff ~ 30, 000 K) may also depend on the assumed mass 
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of the models (currently 0.48 A/©), but this remains to be 
investigated. We find also that stability does depend on 
the assumed surface gravity. Unstable models have log 
g^j 5.7, with a weak dependence on metallicity, but this 
remains also to be investigated in more detail. 

We end this section by presenting, in Table ||, typical 
examples of unstable modes. We list there the periods and 
e-folding times of the unstable modes found in 2 mod- 
els, one computed with Z = 0.06, and the other with Z = 
0.10. Both have log g = 5.2 and log T e g = 4.42. For given 
modes, the periods are significantly larger than in Table ^ 
because of the lower surface gravity considered here (al- 
though for the g-modes the differences are partly due also 
to our use of envelope models instead of full stellar mod- 
els). As expected, the driving is clearly more efficient in 
the high-metallicity model; a wider band of modes are ex- 
cited in this particular model as compared to its Z = 0.06 
counterpart. In both cases, the excited modes have low 
values of the radial order k, and radial as well as nonra- 
dial pulsations are predicted. Since the e-folding times are 
generally significantly shorter than the typical evolution- 
ary time of a hot B subdwarf (~ 10 s yrs), these linear 
instabilities would normally develop into observable am- 
plitudes. 

3. DISCUSSION AND CONCLUSION 

The results presented in this paper pave the way to the 
exciting possibility of being able to use asteroseismology 
as a probe of the internal constitution of subdwarf B stars. 
A subgroup of these stars may indeed constitute a new 
class of pulsating stars. Of course, this possibility hinges 
here on the question of metal enrichment in the driving 
region. We emphasize the fact that this enrichment needs 
only occur in the driving region itself and not necessarily 
everywhere in the envelope (as we assumed in the rela- 
tively crude envelope models used in the present study). 
That heavy elements can show local enrichments in cer- 
tain regions of a hot B subdwarf envelope, and local deple- 
tions in others, is not only plausible but is also expected. 
Indeed, subdwarf B stars all show peculiar surface abun- 
dances whose study constitutes a most active and inter- 
esting subfield of stellar physics in its own right. It is 
currently believed that the abundance anomalies observed 
in the atmospheres of these stars result from the compe- 
tition between gravitati onal settling, radiative levi tation, 
and weak stellar winds ( Michaud et al. 1985 , 1989| ; Lam- 



computations of radiative levitation on metals in models 
of hot B subdwarfs using the same tools considered pre- 
viously in the white dwarf context by Chayer, Fontaine, 
& Wesemael 1995 and |Chayer et al. 1995| . Of high rele- 
vance here, they found that metals, and in particular iron 
(a major contributor to the opacity), do levitate in the en- 
velopes of these stars. Local overabundances of more than 
an order of magnitude can be built through this process. In 
the case of iron, the maximum overabundance occurs near 
log q ~ —9, —10. This result greatly enhances the plau- 
sibility that driving occurs as a result of a local enhance- 
ment of the metal content^ Of course, a detailed study 
of the driving process in hot B subdwarfs will ultimately 
require the use of more detailed envelope structures tak- 
ing into account diffusion processes. In the meantime, we 
feel confident enough to risk the prediction that some sub- 
dwarf B stars should show luminosity variations resulting 
from pulsational instabilities. We believe that the mate- 
rial presented in this paper warrants further theoretical in- 
vestigations of the pulsation properties of these stars and, 
above all, searches for luminosity variations. We have un- 
dertaken both. 

This work was supported in part by the NSREC of 
Canada and by the fund FCAR (Quebec). We thank Matt 
Wood for providing the original versions of the conductive 
opacity subroutines. B. D. acknowledges partial support 
from NASA RTOP 188-41-51-03. 



ontagne et al. 1985, |1987fc |Heber 1987| , |1991| ; Bergeron 
et al. 1988; Chayer et al. 1996 ). In the most recent 
investigation, Chayer et al. 1996 have carried out new 



3 A similar phenemenon may be responsible for mode excitation in rapidly 
oscillating Ap stars. 
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Table 2 

Unstable Modes in Typical Subdwarf B Star Models 
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Fig. 1. — (a) Run of the thermal timescale (expressed in seconds) vs fractional mass depth in a typical subdwarf B model 
(solid curve). The locations of the He burning core, the residual H burning shell, and the Hell-Helll convection zone are 
indicated. As a comparison, the values of the periods of a few low radial order modes with I = 2 are shown. The temperature 
distribution (long-dashed curve) is also illustrated, (b) Run of the Rosseland opacity (dashed curve) and of the integrand of 
the work integral for the / mode with I = 2 (solid curve). The driving region (positive values of dWj dr) is clearly associated 
with an opacity bump, itself caused by heavy element ionization. 
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